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The specific shape of the squark, slepton and gaugino mass spectra, if measured with sufficient 
accuracy, can provide invaluable information not only about the dynamics underpinning their origin 
at some very high scale such as the unification scale Mq, but also about the intermediate scale 
physics encountered throughout their RGE evolution down to the energy scale accessible for the 
LHC. In this work, we study general features of the TeV scale soft SUSY breaking parameters 
stemming from a generic mSugra configuration within certain classes of SUSY 5*0(10) GUTs with 
different intermediate symmetries below Mq- We show that particular combinations of soft masses 
show characteristic deviations from the mSugra limit in different models and thus, potentially, allow 
to distinguish between these, even if the new intermediate scales are outside the energy range probed 
at accelerators. We also compare our results to those obtained for the three minimal seesaw models 
with mSugra boundary conditions and discuss the main differences between those and our 5*0(10) 
based models. 
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I. INTRODUCTION 

All proposed supersymmetry (SUSY) breaking 
schemes have to introduce some high energy scale, where 
soft terms are generated. This scale could be as high 
as the scale of Grand unification (GUT), or even the 
Planck scale in gravity mediated schemes [ll 12], or as 
low as a 100 TeV, for example in gauge mediated SUSY 
breaking (GMSB) [i|. SUSY particle masses at the 
electro-weak (TeV) scale then have to be calculated 
from the fundamental parameters of the models using 
renormalization group equations (RGEs). Although 
those fundamental parameters are a priori unknown, at 
least in minimal schemes there exist certain sum rules for 
SUSY particle masses, which allow to test the different 
SUSY-breaking mechanisms, as has been shown for 
the example of minimal Supergravity (mSugra) already 
some time ago [j] . Based on the detailed studies of the 
capabilities of the LHC and ILC experiments to measure 
SUSY particle masses [5|47|, the accuracy with which 
different SUSY-breaking schemes can be tested has since 
then been calculated by a number of authors, for a few 
examples see J8l4l2|. 



the Minimal Supersymmetric extension of the Standard 
Model (MSSM). Evolution under RGEs is, however, sen- 
sitive to the particle content of the theory. Thus, in prin- 
ciple, any superfield beyond the MSSM (with mass be- 
low the SUSY-breaking scale) will leave its imprint on 
the soft parameters. The specific shape of the squark, 
slepton and gaugino mass spectra, if measured with suf- 
ficient precision, can therefore provide invaluable infor- 
mation not only about the dynamics underpinning their 
origin, but also about physics at intermediate scales. In 
this paper, we study soft SUSY-breaking masses within 
certain classes of SUSY 5*0(10) theories with different 
intermediate symmetries below the GUT scale Mq. Our 
main motivation to study these models comes from the 
observed neutrino masses 13l-lla| and the possibility that 



However, most of these studies concentrated on mod- 
els in which the particle spectrum between the electro- 
weak and the SUSY-breaking scale was exactly that of 



supersymmetry might be discovered soon at the LHC. 

In the MSSM, if SUSY particles have TeV-scale 
masses^, the gauge couplings unify (nearly) perfectly at 
around Mq ^ 2x 10^^ GeV. Adding new particles beyond 
the MSSM spectrum can easily spoil this attractive fea- 
ture and, thus, the requirement of gauge coupling unifica- 
tion (GCU) imposes a severe constraint on SUSY model 
building. However, neutrino oscillation experiments 13- 
ISl . llSJ l have shown that at least two neutrino masses are 
non-zero [l9|, |20[ and at least one neutrino must have a 
mass TOAtm ^ 0.05 eV. If neutrinos are Majorana par- 
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^ Strictly speaking, within SUSY unification requires only gaugi- 
nos to be light [iM. llTll . but not necessarily sferniions. We will 
not entertain this possibility. 



tides, this value indicates that the scale of Icpton num- 
ber violation (LNV), Alnv; can not be larger (but could 
potentially be much smaller, see below) than roughly 
Alnv ~ 10^^ GcV. This value is significantly below Mq. 

In the minimal SUSY SU{5) model neutrinos are mass- 
less, just as in the MSSM and for the same reasons. How- 
ever, it is fairly straightforward to extend minimal SU{5) 
to include a seesaw mechanism which allows to explain 
the observed smallness of the neutrino masses. It is well 
known that, at the renormalizable level, there are ex- 
actly three ways 2l| to do so: (i) Add (at least two) 
gauge singlet superfields, i.e. "right-handed neutrinos", 
this is now usually called type-I seesaw [22 4241 : fii) add a 
scalar triplet with F = 2 (type-II seesaw) 25|, |26[ ; or (iii) 
add (two or more) fermionic triplets with Y = (type- 
Ill seesaw) 27[. For the latter two cases, the successful 
unification of the MSSM can only be maintained if these 
heavy fields enter in complete SU{5) multiplets. Thus, 
within SUSY models, GCU requires the type-II seesaw 
to be realized by adding a pair of Higgs 15-plets, while 
a type-Ill seesaw can be generated with (at least two) 
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copies of 24-plets in the matter sector 

If we require Yukawas to be perturbative anywhere be- 
tween the seesaw scale (i.e., Alnv) and the GUT scale 
Mq, all three types of seesaw require Alnv to be below 
10^^ GeV. If we ask in addition that all gauge couplings 
remain perturbative, lower limits on SUSY type-II seesaw 
of the order of Alnv > 10^ GeV at 1-loop (or Alnv > 10^ 
GeV at 2-loop) [28| result. For type-Ill seesaw, perturba- 
tivity puts lower bounds on the seesaw scale of the order 
of Alnv > 10^^ GeV for three copies of 24 of SU{5) and 
around Alnv ^ 10^ GeV for two^ copies of 24. Since 
type-I seesaw adds only Standard Model (SM) singlets, 
there is no lower limit on its scale from perturbativity 

Models based on S'O(IO) are different from SU{5) 
in that they automatically contain the necessary ingredi- 
ents to generate non-zero neutrino masses: (i) The spino- 
rial 16 of SO{10) contains a complete SM family plus a 
gauge singlet, i.e. a right-handed neutrino. In addition, 
(ii) U{1)b-l is a subgroup of SO{10). If the U{1)b-l 
is broken by SU{2)ji triplets with B — L ^ 2, a see- 
saw mechanism of either type-I and/or type-II results 
automatically [2J1. Alternatively, breaking U{1)b~l by 



^ With only one copy of 24, the seesaw scale could be lowered as far 
as the electro-weak scale. However, with only one 24, neutrino 
data can not be explained unless non-renormalizablo operators 
are added to the model [2911. 



SU{2)fi doublets can give different realizations of the so- 
called inverse [3]| and linear [32| seesaw schemes. 

The SO{10) gauge symmetry can be broken to the SM 
gauge group in a variety of ways [33[ . Since our main mo- 
tivation is neutrino masses, all breaking chains of interest 
to us contain a left-right symmetry (LR) at some stage. 
SUSY LR models which use triplets to break SU{2)ji, 
whether using only B ~ L = 2 triplets 1341 |35J or both 
B — L = 2 and B — L = triplets 36J, |37| , all require 
that the scale of SU{2)ii breaking [vn] is close to the 
GUT scale, typically vr > IQi^ GeV from GCU. How- 
ever, also in non-minimal versions of triplet-based mod- 
els, one can not lower vr arbitrarily, since one encounters 
either problems with proton decay or with perturbativity 
[38|. Allowing for either (a) sizeable GUT scale thresh- 
old corrections, (b) non-renormalizable operators or (c) 
adding some carefully chosen particles the authors of [39 1 
find a lower limit on vr of the order of vr ^ 10^ GeV. 

However, the situation is different in models with dou- 
blets. It was shown in [40| that if one breaks SU{2)fi x 
U{1)b-l by means of an B — L = triplet to C/(l)fl x 
U{1)b-l and, subsequently, the U{1)b x U{1)b~l sym- 
metry gets broken down to U{1)y by the K-neutral com- 
ponents of SU{2)ii doublets, it is possible to construct 
models in which the scale of U(])h x U{1)b-l breaking, 
vbl, can be as low as TeV. In [39| it was demonstrated 
that even the full SU{2)fi can be brought down to the 
electro-weak scale, if only doublets are used in the sym- 
metry breaking and if some additional particles are also 
light. An especially simple variant for a low vr scale was 
discussed in 4l|. Here, GCU is maintained for a TeV 
scale SU{2)ji with only two requirements: (a) The num- 
bers of light left and right doublets have to be different 
and (b) a (pair of) light coloured SU{2)l singlets needs 
to be added to the spectrum. 

Obviously, all models with additional gauge groups 
lead to a potentially very rich phenomenology at the 
LHC. Current limits on new Z' (and Wr) gauge bosons 
are very roughly of the order of mz' ^ (5 — 6)/.g TeV 



{mwn ^ 1 TeV) 42, |43[ , with exact numbers depending 
on the couplings, so there is ample room for discovery. 
One expects that for y/s = 14 TeV at the LHC limits 
for Z' bosons will improve to at least 3 TeV 
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Wr should be discovered at the LHC up to masses of 
the order of 4 - 4.6 TeV 



45 



46|, depending on luminos- 



ity. However, even if the new gauge bosons predicted in 
the models 391441 1 are out of reach for the LHC, sparti- 
cle mass spectra will contain indirect hints for these new 



scales due to changes in the RGEs, as discussed above. 
This observations is in fact the main motivation for the 
calculations presented in this paper. 

Within the niSugra framework, one can define certain 
combinations of soft parameters, which are independent 
of the high scale input parameters at leading order. We 
will call such combinations "RGE invariants" . In l47| it 
was pointed out that these invariants show a character- 
istic deviation from their mSugra expectations, if either 
a type-II or a typc-III seesaw mediators are added to the 
MSSM spectrum. Here, we will study these invariants in 
different 50(10) based models. We wih construct vari- 

and will also con- 
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ants of the models proposed in [40|, 
sider a completely new model, in which vr can be brought 
down to the electro- weak scale with the help of an inter- 
mediate Pati-Salam scale 481. We will show how the 



RGE invariants calculated within these models depart 
from their mSugra values, how they differ from model to 
model and, importantly, also differ from the expectations 
for the minimal type-II and type-Ill seesaws. The invari- 
ants are therefore good indicators to distinguish between 
different GUT-based SUSY models. 

Two comments might be in order at this point. 
First, our calculations rely on the assumption of strict 
mSugra boundary conditions. In principle, invariants can 
be calculated also in more complicated SUSY-breaking 
schemes, if the SUSY-breaking scale is larger than the 
mass scale of the new states. However, with a total of 
only four invariants (per generation) only SUSY-breaking 
schemes with very few additional parameters will lead to 
non-trivial consistency tests. Furthermore, while the in- 
variants are certainly useful model discriminators, it has 
been shown that quantitatively important 2-loop correc- 
tions exist for both, the type-II [49| and the type-Ill see- 
saw 28[ . A quantitative determination of the new inter- 



mediate scales will therefore most likely rely on a detailed 
numerical x^-analysis of measured SUSY spectra [50| , us- 
ing invariants only as guidance for which models might 
be interesting for further scrutiny. 

The rest of the paper is organized as follows. In the 
next section, we shall specify four basic 5*0(10) mod- 
els of interest paying particular attention to their po- 
tential compatibility with the SM flavour structure. In 
section HI, we briefly comment on the evolution of the 
soft masses in mSugra models and, for completeness, re- 
call the deflnitions of the RGE invariants, following es- 
sentially the discussion in [49|. The results of a simple 
analysis of their sensitivity to the intermediate scales in 



the four scenarios considered here are given in Section 
IV. Finally we close with a short discussion and outlook. 
Some technical details of the RGEs in models with more 
than a single abelian gauge factors are deferred to an 
Appendix. 

II. SPECIFIC SUSY SO(IO) GUT MODELS 

Let us begin with a detailed speciflcation of the four 
basic SUSY 50(10) GUTs which shall be the studied 
in Sect. IV. Though all of them, by construction, accom- 
modate the low-energy measured values of the gauge cou- 
plings, they will in general yield vastly different MSSM 
soft spectra whose shapes would strongly depend on the 
character of the intermediate symmetries and the scales 
of their spontaneous breakdown. 

A. General remarks 

In all cases, we demand that the models should be re- 
alistic in several basic aspects and potentially interesting 
for our scope, namely: 

• Requirement 1: SUSY 50(10) unification with a 
sliding intermediate scale by which we mean that 
the position of a certain intermediate scale can be 
moved over a large energy range whilst the full com- 
patibility with the electroweak constraints is main- 
tained. This is a basic practical stipulation in order 
to be able to study the scale-dependence of the soft 
leading-log RGE invariants in such GUTs over a 
large range. 

• Requirement 2: Renormahzable 50(10) -^ MSSM 
gauge symmetry breaking - this is namely to have 
a good grip on the intermediate scales and the as- 
sociated thresholds. 

• Requirement 3: Potentially realistic fermionic spec- 
tra - we demand that the effective Yukawa struc- 
ture is rich enough to be able, at least in principle, 
to accommodate the low-energy matter-fermionic 
spectra and mixing. The sliding nature of the 
SU[2)r X U{1)b-l scale, however, typicall y c alls 
for a non-canonical seesaw, such as inverse 3l[ or 
linear J32l | seesaw. 

• Requirement 4: MSSM Higgs doublet structure 
suitable for the implementation of the standard ra- 
diative symmetry breaking and also as a means 



to get unrelated Yukawa couplings for quarks and 
charged leptons. 

As to the Requirement 1 above, we shall be namely in- 
terested in SUSY SO{10) models with a shding SU{2)r 
breaking scale which would be assumed to range from 
as low as several TcV up to essentially the GUT scale. 
From the gauge unification perspective, there are two ba- 
sic strategies to devise such Models. In practice: 

• One can attempt to compensate for the departure 
of the fo-coefficients from their "canonical" MSSM 
values (due to the presence of W^ and the SU{2)]i- 
breaking Higgs multiplets in the desert) by other 
multiplets brought down to the S'C/(2)ij-breaking 
scale, which would inflict further shifts to the h- 
coefficients (namely 173) in order to compensate for 
the genuine low-scale SU{2)ji effects. The main ad- 
vantage of this approach is that 5'C/(2)/j x U{1)b-l 
becomes the only intermediate scale at play, so the 
50(10) gauge symmetry is broken down to the 
SU{3)c X SU{2)l X U{1)y of the MSSM in just 
two steps. The slight complication here is the fact 
that the gauge-coupling unification in such a case 
is not exact, which brings an extra theoretical un- 
certainty into the game. ^ 

• Alternatively, rather than compensating for the 
departure of the 6-coef5cients from their MSSM 
values due to the Wj^ (and the associated Higgs 
multiplets) in the desert, one can take advan- 
tage of this and invoke an extra intermediate scale 
such as for instance SU{4)c x SU{2)l x SU{2)r 
of Pati and Salam and let it conspire with the 
SU{2)r X U{1)b-l so that the gauge unification 
is maintained. Though this can be somewhat more 
elaborate in practice, the clear advantage of such 
a scenario is that one can always devise an exact 
gauge coupling unification by a proper adjustment 
of the Pati-Salam scale. 

In both cases, because we have quite a lot of beyond- 
MSSM dynamics in the desert, we expect significant ef- 
fects of the relevant intermediate scale(s) on the shape of 
the MSSM squark and slepton spectra. 

The first strategy above, especially in combination 
with the other requirements, is rather restrictive. In- 
deed, it imposes strict conditions on the 6-coefficients in 



specific models which should essentially match those of 
the MSSM up to a uniform shift. Nevertheless, a variety 
of potentially realistic models can still be devised and, in 
particular, the behaviour of the RGE invariants in this 
class of theories can be strongly model-specific. We shall 
demonstrate this on a couple of scenarios of this kind 
derived from [4l|, c.f.. Model I and Model II in section 

EH 

The sensitivity to the intermediate-scales dynamics 
should be even more pronounced in the latter class of 
scenarios with more than a single such scale at play. 
This is namely due to the fact that the extra fields 
in the desert associated to a higher intermediate sym- 
metry (e.g., Pati-Salam) tend to affect the soft spectra 
stronger than in the former case with an intermediate 
SU{3)c X SUi2)L X SU{2)r x U{1)r only. This feature 
is going to be clearly visible in the specific model of this 
kind, c.f.. Model III in section Hi Bl 

However, a strong dependence of the invariants on the 
sliding scale should not be viewed as a generic feature of 
the SUSY 50(10) GUTs. Indeed, there are simple sce- 
narios in which the sliding intermediate scale does leave 
almost no imprints in the soft spectrum. We shall demon- 
strate this on a specific model with a sliding intermediate 
U{1)b-l scale (and a fixed SU{2)r scale at around lO^'^ 
GeV ensuring a proper gauge unification) of the kind 
given in [40|, c.f.. Model IV in the section Hi Bl Here, the 
GUT-scale pattern of the RGE invariants is (almost) not 
changed by the running, leaving no good handle on the 
intermediate scale in the SUSY spectrum. 



B. SUSY 50(10) models with a sliding SU{2)r scale 

1. Models I and II: single sliding intermediate scale 

First, we shall introduce two variants of the model ad- 
vocated in 



4l| which supply the original setting with a 



This feature is already present at the MSSM level. 



few extra ingredients in order to make it potentially re- 
alistic, c.f.. Sect. Ill Al 

a. Model I: The field content relevant to the running 
in this Model is specified in TABLE HI 

The original 50(10) gauge symmetry is broken down 
to the MSSM in two steps via an intermediate SU{3)c x 
SU{2)l X SU{2)r X U{1)b-l symmetry stage. The first 
step, i.e., the 50(10) -^ SU{3)c x SU{2)l x SU{2)r x 
U{1)b-l breaking, is triggered by an interplay between 
the VEVs of 45 and 210 of 50(10). Subsequently, the 



Field 


Multiplicity 


3c2l2r1b^l 


50(10) origin 


Q 


3 


(3,2,1,+i) 


16 


Q^ 


3 


(3,1,2,-i) 


16 


L 


3 


(1,2,1,-1) 


16 


L" 


3 


(1,1,2,+1) 


16 


S 


3 


(1,1,1,0) 


1 


Sd, 5d 


1 


(3,1,1,-1), (3,1, 1,+f) 


10 


$ 


1 


(1,2,2,0) 


10, 120 


X, X 


1 


(1,2,1,±1) 


T6, 16 


x^r 


3 


(l,l,2,Tl) 


T6, 16 



TABLE I: The relevant part of the field content of Model I 
with a sliding S'l7(2)fl;-breaking scale discussed in section 
111 B 1 al In the third column the relevant fields are character- 
ized by their SU{3)c x SU{2)l x SU{2)r x U{1)b-l quantum 
numbers while their S'O(IO) origin is specified in the fourth 
column. 



SU{2)f{ X U{1)b-l gEi'Ugc symmetry is broken down to 
the U{1)y of hypereharge by the VEVs of the x'^ © X'^ 
pair which can emerge for instance from its cubic in- 
teraction with a fuU singlet p^. Note that, at the one- 
loop level, such a neutral p field can be put essentially 
anywhere between Mqut and Msusy without any im- 
pact on the gauge unification. The 6i-coefficients at the 
SU(2,)c X SU{2)l X SU(2)rx U{1)b-l level read b^ = -2, 
bL = 2, bft = 4: and 6^_£, — 13 where the last num- 
ber corresponds to the canonically normalized B — L 
charge, which is obtained from the "physical" one (based 
on Bq = +-5 and L^ ~ +1) by means of the formula 

(B - Lf = J^{B ~ L)P. Note that these coefficients 
happen to be entirely identical to the setting advocated 
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and, thus, the leading-log evolution of the soft 
SUSY-breaking parameters in that Model is covered by 
our analysis of Model I, see Sect. IIV Al The scale of 
the SU{2)ji X U{1)b-l breaking is not determined be- 
cause it drops from the formula for the unification scale 
(owing namely to the hypercharge-matching condition 
Uy^ = 5 Q^i?^ + Z'^b^-l) ^^"^ affects only the value of the 
GUT-scale gauge coupling ac which, however, is subject 
of much weaker constraints. 

This, on one hand, makes the gauge coupling unifica- 



* It is perhaps worth mentioning that, for a very low scale of 
the SU(2)ji X U{1)b-l breaking, the relevant VEV can be de- 
vised even without an extra singlet because, then, the inter- 
play between a "RH /^-torm" for x'^ ffi X'^ aid the relevant soft 
mass should be sufficient, in complete analogy with the SU{2)l- 
doublet sector in the MSSM. 




FIG. 1: Gauge coupling unification in Model I in two limits 
corresponding to different positions of the sliding SU{2)r x 
U{1)b-l breaking scale Vr. In solid fines, we depict the RGB 
behaviour of the gauge couplings for Vr in the vicinity of 
the electroweak scale Vr ~ lO^GeV while the dashed lines 
correspond to Vr ~ 10^*GeV. The position of the intersection 
region shifts slightly up with rising Vr but the corresponding 
scale remains intact. 



tion in Model I qualitatively similar to the MSSM case, 
see FIG. [TJ On the other hand, it is well known that in 
the MSSM the one-loop gauge unification is incompatible 
with the latest extractions of as unless the soft SUSY- 
breaking scale is pushed well below 1 TeV. This few per- 
cent mismatch is expected to be accounted for by GUT- 
scale thresholds whose detailed analysis is, however, be- 
yond the scope of this work. Thus, in what follows, 
we shall simply parametrize our ignorance of the shape 
of the GUT spectrum by considering unification regions 
from where the SU{3)c x SU{2)l x SU{2)r x U{1)b-l 
gauge couplings can emanate rather than unique unifi- 
cation points, c.f., FIG. [5] and discussion in Sect. IIVAI 
Though this approach is oversimplified in several aspects, 
it admits to estimate the magnitude of the theoretical 
error associated to the lack of exact gauge-coupling uni- 
fication in this model. 

Concerning the effective flavour structure of the model, 
there are two aspects worth some discussion here, namely, 
the structure of the effective MSSM Higgs doublet pair 
$ and the possibility to accommodate the SM quark and 
lepton masses and mixing (requirements 3 and 4 formu- 
lated at the beginning of this Section). First, the effective 
L-R bidoublet (1, 2, 2, 0) corresponds to a massless com- 
bination of the (1, 2, 2) and (15, 2, 2) Pati-Salam compo- 
nents of 10 and 120 of 50(10), respectively, which can 
mix at the GUT-lcvel due to the PS-breaking VEV in an 
S'O(10)-brcaking multiplet such as 45 and/or 210. Usu- 



Field 


Multiplicity 


3c2l2r1b^l 


50(10) origin 


Q 


3 


(3,2,1,+i) 


16 


Q^ 


3 


(3,1,2,-i) 


16 


L 


3 


(1,2,1,-1) 


16 


L" 


3 


(1,1,2,+1) 


16 


S 


3 


(1,1,1,0) 


1 


Su, Su 


1 


(3, 1,1,+|), (3,1,1,-1) 


45 


$ 


1 


(1,2,2,0) 


10, 120 


X, X 


1 


(1,2,1,±1) 


T6, 16 


x^r 


2 


(l,l,2,Tl) 


T6, 16 




TABLE II: The same as in TABLE U for Model II defined in 
Sect. HTB 1 bl The main variation with respect to Model I is 
the B — L charge of the vector-like colour triplet pair owing 
to its different 50(10) origin. The extra Su and 5u fields can 
mix with the up-type quarks at the MSSM level which leads 
to a potentially realistic effective flavour structure. In order 
to maintain the MSSM-like unification pattern, the number 
of the SU{2)r doublets has been reduced, thus making the 
setting slightly more compact than in Model 1. 



FIG. 2: Gauge couphng unification in Model II in two limits 
corresponding to different positions of the sliding SU{2)r x 
U{1)b-l breaking scale Vr. In solid fines, we depict the RGE 
behaviour of the gauge coupfings for Vr in the vicinity of 
the electroweak scale Vr ~ lO^GeV while the dashed lines 
correspond to Vr ~ lO^^GeV. The position of the intersection 
region shifts slightly up with rising Vr but the corresponding 
scale remains intact. 



ally, the role of the extra Higgs such as 120 and/or 126 in 
the Yukawa sector is namely to provide Clebsch-Gordan 
coefficients that would break the degeneracy of the effec- 
tive Yukawa couplings among up and down quarks and 
charged leptons. However, an extra 120 alone is still not 
enough as it does not yield enough freedom to accommo- 
date the SM data [55|. Actually, the issue becomes even 



worse if the MSSM-level mass matrices for the two hy- 
percharge components of $ arc virtually identical, as one 
can expect for a single bidoublet at play in the low-scale 
S'C/(2)ij-breaking regime. Both issues are potentially re- 
solved due to the extra vector-like down-type quark pair 
Sd®Sd and an additional S'[/(2)i-doublet Higgs pair x®X 
(which, simultaneously, ensure the right &-coefficients for 
the running), c.f., TABLE HI In this case, the down- 
type quark mass matrix is extended^ to 4 x 4 which, 
together with the extra freedom in the MSSM doublet 
sector, should be enough to avoid the Grimus-Kuhbock- 
Lavoura (GKL) no-go 



551. Let us also mention that the 



VEV of XL gives rise to the LS entry in the neutrino mass 
matrix gen erating the linear seesaw mechanism and, un- 
like in [40|, it is not naturally suppressed in this case 
because x © x resides well below the GUT scale. Thus, 
one has to assume a small LSx Yukawa coupling. 



^ For an explicit 50(10) realisation of this mechanism see e.g. [56 
and references therein. 



b. Model II: The relevant 6i-coefficients at the 
S't/(3)cxS'C/(2)Lx5t/(2)KxC/(l)B_L level read 63 = -2, 
6l = 2, bn = 3 and 5f^^ = 29/2. Indeed, these num- 
bers differ from Model I only in the SU{2)r x U{1)b-l 
sector and the variations in the relevant fo-coefficients 
obey Abu + I^^b-l = so the 6-coefficient associ- 
ated to the "effective" MSSM hypercharge is the same 
as in Model I. Therefore, apart from the difference in 
the specific slopes of the SU{2)ii x U{1)b-l curves the 
qualitative picture of the gauge coupling unification in 
Model II, c.f. FIG. [21 is very similar to that observed 
in Model I. Nevertheless, as we shall see in Sect. IIV Al 
even such a slight change in the gauge-coupling behaviour 
at the SU{3)c x S'C/(2)l x SU{2)ii x U{1)b~l level is 
enough to generate a significant difference between the 
Model-I and Model-II soft invariants, especially if the 
SU{3)c X SU{2)l X SU{2)r x U{1)b-l running is long. 
However, if the SU{2)ji x U{1)b-l gauge symmetry hap- 
pens to be broken close to the GUT scale, the two mod- 
els will be indistinguishable from the soft-sector point of 
view. 

Concerning the flavour structure of Model II, it is in- 
deed very similar to that of Model L with the main dif- 
ference that here the GKL no-go [55[ is overcome by a 
4x4 extension of the up-type quark mass matrix. More- 
over, since it is the VEV of x rather than that of x that 
enters the extended up-type quark mass matrix, (x) can 
be made much smaller than (x) which also relieves the 
need for the small LSx Yukawa in the neutrino sector. 



Field 


Mult. 


3c2l2r1b-l 


Pati-Salam 


50(10) 


Q 


3 


(3,2,1,+i) 


(4,2,1) 


16 


Q^ 


3 


(3,1,2,-i) 


(4,1,2) 


16 


L 


3 


(1,2,1,-1) 


(4,2,1) 


16 


L" 


3 


(1,1,2,+1) 


(4,1,2) 


16 


E'^ 


3 


(1,1,3,0) 


(1,1,3) 


45 


Sd, Sd 


1 


(3,l,l,Ti) 


(6,1,1) 


10 


$ 


2 


(1,2,2,0) 


(1,2,2) 


10 


n 


1 


(1,1,3,0) 


(1,1,3) 


45 


x,x 


1 


(1,2,1,±1) 


(4,2, 1), (4, 2,1) 


16, 16 


x^r 


1 


(l,l,2,Tl) 


(4,1,2),(4,1,2) 


16, 16 


* 


1 


absent 


(15,1,1) 


45 



TABLE III: The eflective field contents of Model III in the 
two intermediate symmetry stages. 

Given also the reduced number of the SU{2)ii doublets, 
Model II constitutes a somewhat more compact alterna- 
tive to Model I. 




FIG. 3: Running in the Model III variant of the low-LR scale 
SUSY SO(IO). Here the SO(IO) gauge symmetry is broken 
first into a Pati-Salam intermediate stage stretching from the 
unification point down to the relevant energy scale Vps (in the 
middle) and, subsequently, to the L-R symmetry stage. The 
value of Vps is correlated to the position of the L-R breaking 
scale Vr which can again slide from as low as few TeV up to 
roughly 10'*GeV, c.f., FIG. El 



2. Model III: sliding SU{2)r and Pati-Salam scales 

The third model of our interest belongs to the second 
category of the simple classification given in Sect. Ill Al In 
particular, the sliding nature of the SU{2)b, x U{1)b-l 
scale is achieved via an interplay with another interme- 
diate scale, namely, the Pati-Salam SU{4:)c' x SU{2)l x 
SU{2)f{, rather than a delicate adjustment a la Model I 
or Model II owing to very specific field contents. Thus, 
the initial SO{10) gauge symmetry is broken down to the 
MSSM in three steps. The field content relevant to the 
two intermediate-symmetry stages is given in TABLE lTlll 
In more detail, the initial 50(10) -^ SU{A)c x SU{2)l x 
SU{2)r breaking is triggered by the GUT-scale VEV 
of 54 of 50(10). The subsequent SU{A)c x SU{2)l x 
SU{2)r -^ SU{3)c X SU{2)l x SU(2)r x U{1)b-l break- 
ing is due to the VEV of the ^ emerging again from its 
interplay with an extra singlet. Finally, the SU{2)r x 
U{1)b-l symmetry is broken down to the MSSM by 
means of the VEVs of x'^ © x'^ which are connected by 
the B - L-neutral S'C/(2)fl-triplet fl. At the Pati-Salam 
stage, the 6i-coefficients read 64 ~ 3, b^ = 6, bji = 14 
while at the SU{3)c x SU{2)l x SU{2)r x U{1)b-l level 
they are 63 = -2, bL = 3, bR = 11 and 6^^f ^ = 10. 

In this model, both the position of the GUT scale as 
well as the value of ac depend on both intermediate 
scales. However, unlike in Models I and II, here the gauge 
unification can always be made exact, c.f., FIG.|31 even at 



the one-loop level, and, thus, there is no extra theoretical 
uncertainty other than the error in the electroweak-scale 
Og to be taken into account. 

The flavour structure of this model relies on the pres- 
ence of three extra copies of SU{2)r triplet S'^ which in 
the neutrino sector play a role similar to that of S in 
Models I and II. In particular, they expand the 6x6 
neutrino mass matrix to 9 x 9 where, e.g., the L'^E'^ sec- 
tor comes from the contraction with the VEV of x'^, but 
without any entry generated at the LT,'^ "linear seesaw" 
position. Thus, there is no need for an extra fine-tuning 
in the seesaw formula in Model III. Moreover, the charged 
components of T,'^ can mix with the charged leptons and, 
hence, provide the welcome departure from the down- 
type quarks even if the MSSM doublets span over lO's of 
50(10) only. Indeed, the relevant 6x6 charged-lepton 
mass matrix looks schematically like 



Mi ex 



Yv 



10 







(r> Ms^ + (f^> 



(1) 



where the row and column bases are {L, S"^^} and 
{L'^,I]'^+}, respectively, and fi-^t: is the associated sin- 
glet mass parameter. Note also that the VEV of fl is 
antisymmetric in the generation space and, thus, does 
not contribute to the neutrino Majorana mass matrix. 
Finally, the two MSSM Higgs doublets arc different be- 
cause the underlying bi-doublets contract through 51 and, 
therefore, the effective up-type quark Yukawa coupling 
differs from the down-type one even without the need to 



resort to the mixing with the vector-hke Sd © Sd pair. 



C. SUSY S'O(IO) models with a sliding U{1)r scale 

All the models discussed so far featured an intermedi- 
ate SU{2)r X U{1)b-l symmetry which, at a certain 
scale, was broken directly down to the U{\)y of the 
MSSM hypercharge. The full SU{2)r, however, is not 
the minimal option to realize a gauge symmetry acting 
in the RH sector of the matter spectrum in a way com- 
patible with the MSSM quantum numbers. Indeed, the 
hypercharge sum-rule Y = T^ + {B — L)/2 trivially holds 
even if one sticks to the U{l)ii subgroup of the original 
SU{2)ii generated by T^ alone. 

On the other hand, within 50(10) broken down to an 
intermediate SU{3)c x SU{2)l x U{1)r, x U{1)b-l stage, 
only Z' and the associated U{1)r x U{1)b-l — > U{1)y- 
breaking Higgs fields remain light (at least in minimally 
fine-tuned scenarios) and, thus, the intermediate-scale 
dynamics is generally much simpler than in the models 
based on the full SU{2)u x U{1)b-l- In view of that, one 
can expect that also the intermediate-scale dependence of 
the soft RGE invariants will be much milder than in the 
former case. Moreover, with more than a single abelian 
gauge factor at play, there is a new class of effects as- 
sociated with the so called kinetic mixing between the 
associated gauge fields. Both these aspects make this 
class of models worth further scrutiny. 



General remarks 
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Mult. 


■ic2LlRlB-L 


3c2l2r1b-l 


SO(IO) 


Q 


3 


(3,2,0,-fi) 


(3,2,1,-fi) 


16 


Q^ 


3 


(3,l,±i-|) 


(3,1,2,-i) 


16 


L 


3 


(1,2,0,-1) 


(1,2,1,-1) 


16 


L" 


3 


(l,l,±i+l) 


(1,1,2,+1) 


16 


S 


3 


(1,1,0,0) 


(1,1,1,0) 


1 


$ 


2 


(l,2,±i,0) 


(1,2,2,0) 


10 


Vt 


1 


absent 


(1,1,3,0) 


45 


x,x 


1 


absent 


(1,2,1,±1) 


T6, 16 


x\T 


1 


(i,i,±iTi) 


(l,l,2,Tl) 


T6, 16 



TABLE IV: The eflective field contents of Model IV relevant 
to the two intermediate symmetry stages. 

at the one-loop level. Thus, the specific position of the 
U{1)i{XU{1)b-l -^ U{1)y breaking scale is, in this case, 
irrelevant for the one-loop gauge running. 

This, however, is not the case for the leading-log soft 
RGE invariants of our interest. In particular, unlike By, 
both the C/(1)_R x U{1)b-l gauge bosons enter the renor- 
malized propagators of squarks and sleptons 





$ / G 



(2) 
and one can expect a residual dependence of the invari- 
ants on the U{l)ii x ?7(l)B_i-breaking scale. Neverthe- 
less, as we shall demonstrate in a particular realization 
of this simple scheme, such effects should be much milder 
than those in the scenarios with the full gauged SU{2)ji 
symmetry. 



Remarkably, the simplicity of the minimally fine-tuned 
U{l)ji X U{1)b-l ^- U{1)y scenarios automatically im- 
plies the scale of this spontaneous symmetry breakdown 
is a sliding one. Indeed, minimal fine-tuning implies that 
the spectrum of the Model in the unbroken phase con- 
sists of that of the MSSM plus Z' plus an MSSM-neutral 
Higgs responsible for the relevant symmetry breaking. 
Since the gauge field associated to the hypercharge (By) 
does not feel any effect of either Z' nor the hypercharge- 
neutral Higgs^ the "effective" hypercharge gauge cou- 
pling (corresponding to a relevant combination of gy and 
gfj) in this picture runs as if it were in the MSSM, at least 



2. Model IV: U{1)r x U{1)i 



U{1)y breaking 



To put this statement on a firm ground the effects of the kinetic 
mixing must be considered, see, e.g., 1531 . 



Here we consider a variant of the basic SUSY 50(10) 
model advocated in [40[ in which an extended interme- 
diate U{1)r X U{1)b-l stage follows a short SU{2)ji x 
U{\)b-l phase. The field content relevant to the RG 
running in the first two parts of the symmetry-breaking 
chain is given in TABLE |IVl 

In more detail, after the initial 50(10) -^ SU(3)c x 
SU{2)l X 5C/(2)fl X U{1)b-l breaking triggered by es- 
sentially the same mechanism as in Models I and II, the 
subsequent 5C/(2)h x U{1)b-l ^ U{1)b x U{1)b^l re- 
quires a VEV of the O field which, at the level of an 
effective theory, can again emerge from its interplay with 
a LR singlet. The last symmetry-breaking step is then 



60. 
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10'" 
E (GeV) 
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FIG. 4: Running in Model IV with a low B — L scale. Note the 
effects of the U(l) mixing in the running & matching; the low- 
est curve corresponds to the off-diagonals of the (GG"^/47r)~^ 
matrix. 



achieved in a similar manner by the VEVs of the MSSM- 
neutral components of the x'^ © It fields. 

The relevant 6i-coefficicnts at the SU{i)c x SU{2)l x 
5f/(2)fl X Uil)B-L level read 63 = -3, ^l = 2, 6^ - 5 
and 6f^i = 15/2. In the SU{3)c x SU{2)l x U{1)r x 
U{1)b-l stage, however, the effects of the U{1) mixing 
must be taken into account and, thus, the 6-coefficients 
in the t/(l)_R x U{1)b-l sector constitute a matrix of 
anomalous dimensions 7. One has 63 = —3, b^ = 1 and 



\ -1 18 



(3) 



which should be brought into the canonical basis by 
means of a normalization matrix N = diag(l, -^/S/S), 
^can ^ Ar^phys^_ The details of the one-loop RGE evo- 
lution of gauge couplings and soft masses in theories with 
more than a single abelian gauge factor are summarized 
in Appendix [A] The qualitative features of the gauge- 
coupling running in this setting can be seen in FIG. HI 

Concerning the flavour structure of the model, the sit- 
uation is essentially identical to that described in the 
original work [40|. The only exception is a second LR 
bi-doublet retained until the U{1)b. x U{1)b-l breaking 
scale which might be necessary in order to get a poten- 
tially realistic pattern of the effective Yukawa couplings. 
Nevertheless, the salient features of the Model in the soft 
sector should not depend much on the detailed realiza- 
tion of the effective Yukawa pattern. 



III. LEADING-LOG RGE INVARIANTS 

In this section we focus on the calculation of the in- 
variants using mSugra boundary conditions. mSugra is 
defined at the GUT-scale, Mq, by a common gaugino 
mass Ml 12, a common scalar mass ttiq and the trilinear 
coupling Aq, which gets multiplied by the correspond- 
ing Yukawa couplings to obtain the trilinear couplings in 
the soft SUSY breaking Lagrangian ^. In addition, at the 
electro- weak scale, tan/3 = Vu/vd is fixed. Here, as usual, 
Vd and Vu are the vacuum expectation values (vevs) of the 
neutral components oi Hd and 7J„, respectively. Finally, 
the sign of the /i parameter has to be chosen. 

Renormalization group equations for general super- 
symmetric models are known up to 2-loop order |5l| . The 
only case not covered in the otherwise general equations 
given in |5l| are supersymmetric models with more than 
one C/(l) group. With more than a single abelian gauge 
factor, there appears a new class of effects associated 
with the so called kinetic mixing between the associated 
gauge fields. RGEs for this case have been derived very 
recently in 52 1. 

Barring for the moment the effects of C/(l) mixing 
present in Model IV, at the 1-loop level, one can de- 
vise a simple set of analytic equations for the soft terms. 
Gaugino masses scale as gauge couplings do and so the 
requirement of GCU fixes the gaugino masses at the low 
scale 



M.i{msusY) = 



a{M( 



M 



G) 



1/2- 



(4) 



Eq. (HI implies that the relationship of the Mi to M1/2 
is changed in Models I to III, since a{MG) is shifted. 

Neglecting the Yukawa couplings for the soft mass pa- 
rameters of the first two generations of sfcrmions one can 
write 



M 



m 



1/2 



/ 



aiAIc) 



N 

2EE//'«^K)' 

it:, i=i 



(5) 



Here, the sum over "i?j" runs over the different regimes 
in the models under consideration, while the sum over i 
runs over all gauge groups in a given regime, at (va^ ) is to 



'^ It is sometimes argued that this setup should better be called 
CMSSM, since there are even simpler models of supergravity type 
breaking in which Aq is not a free parameter, as for example in 
Polonyi type supergravity 1531 . 1541 . Since we will be concerned 
with only the first two sfermion generations this distinction is 
irrelevant for us. 
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be understood as the gauge coupling of group i evaluated 
at the upper end of regime Rj . In the MSSM one would 
have only to consider one regime, namely from the SUSY 
scale to the GUT scale. In Models I and II we have two 
different regimes, while in Models III and IV there are a 
total of three regimes to consider. 

The different //*" can be written in a compact form as: 

1 



/f - 



ai{vy) 



(6) 



where Vx and Uy, respectively, indicate the value of the 
relevant a at the lower and higher boundaries of the 
regime under consideration. The c^ ' coefBcients given in 
TABLE |V] arc proportional to the values of the quadratic 
Casimir of representation Rf hosting the matter field / 
with respect to the group G in the regime R 

cP = 2CG(i?/). (7) 

They are readily evaluated from the basic formula 

CG{R)d{R)^T2{R)d{G), (8) 

where d(G) is the dimension of the group G, T2{R) the 
Dynkin index of the representation R and d{R) is its 
dimension. Note that the coefficients c^ are differ- 
ent for the different fermions, which leads to a differ- 
ent coefficient in front of M1/2 in eq. ([S]). The bi in 
eq. © arc the one-loop ^-coefficients for the differ- 
ent models defined in the previous section. For com- 
pleteness, the well-known one- loop beta-coefficients for 
the MSSM are (in the traditional SU{5) normalization): 

Eq. (O is valid neglecting C/(l)-mixing effects. The 
extra effects due to the kinetic U{1) mixing relevant in 
Model IV are summarized in the Appendix I3 for a more 
detailed discussion including higher-loop effects, see [52| . 

Individual SUSY masses depend strongly on the initial 
values for toq and Afi/2. However, one can form four 
different combinations, which we choose to be 



LE = (m? 



l)/Ml^ 



(9) 



QE={ml-ml)/Ml 
DL={m%-ml)/Ml 
QU={ml-ml)/Ml 

It is easy to see that, at the leading-log level, mo and 
Ml/2 drop out of the equations for the invariants. Note, 
that one could have equally well normalized to any of 
the other two gaugino masses. The choice of A/i is only 
motivated by the expectation that it will be the gaugino 
parameter measured with the smallest error. 
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L 
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MSSM 




f.MSSM 
^1 
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15 


8 1 
15 30 


f.MSSM 
^2 





3 

2 





1 


f,MSSM 
C3 








8 
3 


8 8 
3 3 


[/(l)fl X U{1)b-l 




J,BL 
'^BL 


3 

4 


3 

4 


1 
12 


1 1 
12 12 


f.BL 





3 

2 





1 


J,BL 


1 

2 





1 

2 


1 


J.BL 








8 
3 


8 8 
3 3 


LR 




J,LR 
^BL 


3 

4 


3 

4 


1 
12 


1 1 
12 12 


f.LR 





3 

2 





1 


J,LR 


3 

2 





3 

2 


1 


^f,LR 








8 
3 


8 8 
3 3 


Pati-Salam 




e^^ 





3 

2 





1 


ci^- 


3 

2 





3 

2 


1 


ci'-' 


15 
4 


15 

4 


15 

4 


15 15 

4 4 



TABLE V: Coefficients c{ for eq. © for difi^erent symmetry 
stages. The MSSM and the LR parts are relevant to all four 
models under consideration; the U{X)r x U{1)b-l and the 
Pati-Salam parts are used solely for Model IV and Model III, 
respectively. 



IV. SLIDING SCALE IMPRINTS IN THE 
LEADING-LOG RGE INVARIANTS 

A. Models I and II with a sliding SU{2)r scale 

The method: As we have already mentioned in 
Sect, ini in Models I and II the sliding nature of the 
SU{2)ii scale makes it impossible to get an exact unifica- 
tion, in full analogy with the MSSM. Since, however, this 
is just about a 2% effect, we shall not attempt to improve 
on this by cither looking for a suitable set of threshold 
corrections or by going beyond the one-loop approxima- 
tion*. Rather than that, we shall just parametrize our 
ignorance of the "true values" of the unification scale po- 
sition and the unified gauge coupling in terms of a pair of 
small "offset" parameters scanning over the area of the 
relevant "non-unification triangle" shown in FIG. ([5]). In 
what follows, we shall to use the error on as{Mz) given 



* Indeed, this would be inconsistent as we are concerned only with 
the leading-log approximation for the softs. 
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FIG. 5: The MSSM-like non-unification triangle in Models I 
and III with vr = 10^* GeV for two different values of the 
unknown soft-SUSY breaking scale (msuSY ~ 1 TeV for the 
upper one and msuSY = 500 GeV for the lower). The upper 
sides of the triangles corresponds to a£^ while the lower- left 
sides depict the "effective" Oy^ defined as |a]^^ + TCtg^^. 
The light blue area surrounding the a^ line represents the 
la uncertainty in aaiMz) as given in [57|]. Both triangles 
move down for lower values of vr, see FIGs. [1] and O 



57|, A{as{Mz)) = 0.002, which docs not take into 



account the latest QCD lattice calculations results. 

The results: In FIGs. |6] and [7] we display the vr- 
dcpendencc of the RGE invariants in Models I and II due 
to the running effects subsumed by Eq. ([5]). The bands 
correspond to the error in the gauge-coupling unification 
inherent to these settings which, at the leading-log level, 
can be taken into account by scanning over the area of the 
relevant non- unification triangle, c.f., FIG.[5l The upper 
(yellow) band refers to the combination QE, the (blue) 
band which at low Wfl partially overlaps with QE rep- 
resents DL, whereas the third (brown) band is QU and, 
finally, the lowest (green) band refers to the LE combina- 
tion. Note that, for practical reasons, the invariants QE 
and DL have been scaled down by a factor of ten. The 
same colour-code is adopted in the other figures in this 
section. 

Several comments are in order here: In general, the 
invariants exhibit a logarithmic dependence on vn. For 
vu close to the MSSM scale (on the left), the QU and 
LE invariants overlap. This is attributed to the enhanced 
gauge symmetry throughout the whole 771susy--^^g range 
which makes m% and m% as well as m? and ni^ behave 

the same, see the LR-stage q -coefficients in TABLE [V] 
In the Vfi — > Mq limit, the mSugra values of the invari- 
ants (modulo the MSSM non-unification) are reproduced. 
Concerning QE and DL, the first thing to notice is that 
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10"' 



10'2 

[GeV] 



FIG. 6; The ii_R-dependence of the leading-log invariants in 
Model I, c.f., Sect. Ill B 1 al The bands represent the error due 
to the non-exact gauge-coupling unification depicted in FIG [5] 
For practical reasons, the numerical values of the invariants 
QE and DL have been scaled down by a factor of ten. 
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FIG. 7; The same as in FIG. E] but for Model II of 
Sect. IIIB Ibl The QU and LE behaviour differs from that 
in FIG. |6] mainly in the 1ow-i;_r regime. 



these invariants tend to increase with vr departing from 
Mq , thus leading to a pattern characteristic to this class 
of models. Moreover, they are more sensitive to the ini- 
tial condition because the colour-effects in their evolution 
do not cancel, thus leading to larger bands. 

Naturally, the main difference between FIG. [7] and 
FIG.[6]is expected in the low-w^ regime where the effects 
due to the slight difference in the Model-I and Model-II 
spectra arc most pronounced and the QU and LE in- 
variants run faster due to a larger ratio of the coupling 
constants in the relevant Eq. ^. 
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FIG. 8: The correlation of the intermediate symmetry- 
breaking scales in Model III (allowed region coloured). The 
contours correspond to the quality of the fit of as{Mz) for 
each choice of the Pati-Salam breaking scale vps an the LR 
breaking scale vr, within la (white area within the coloured 
band), 2a (yellow) and 3a (orange) of the range quoted in 

El. 



B. Model III with sliding SU{2)r and PS scales 

The method: In Model III, the LR and PS interme- 
diate seales can be always adjusted so that one gets an 
exact one-loop unification for vn stretching up to about 
10" GeV, c.f., FIG. E This is technically achieved by 
relating the value of the PS scale to the value of the LR 
scale as 

tps - ^tLB^—{lAtsuSY + 20tz (10) 

+ Tr{18as{tzr^ - 33aL(tz)"^ + 15ay(tz)"^)) 

Here, the tx stand for ln{mx) as usual. Thus, the main 
uncertainty at this level comes from the experimental er- 
ror in as{Mz)- In what follows, we shall vary vr and vps 
along the constant as(A/2 )-error trajectories, namely, 
within ±1(7, corresponding to the boundaries between 
the yellow and white areas within the parameter area 
depicted in FIG. H 

The results: In this case, the intermediate-scale de- 
pendence of the leading-log RGE invariants is yet more 
pronounced than in Models I and II, c.f., FIG[9l As be- 
fore, the numerical values of the invariants QE and DL 
have been conveniently scaled down by a factor of ten. 
For each of the four invariants, the colid curve FIG [9] 
corresponds to as{Mz) fixed at its central value and the 
dashed and dotted lines refer to the — Icr and +lcr tra- 
jectories, respectively. On the high-w^j tail, the different 
curves stop at different energies due to the need to re- 
spect the natural vr < vps hierarchy reflected by the 
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FIG. 9: Intermediate-scale dependence of the RGE invariants 
Model III, see Sect. Ill lT2l For each of the four invariants, the 
solid curve corresponds to as(J^lz) fixed at its central value 
and the dashed and dotted lines refer to the —Icr and +lcr 
trajectories, respectively; c.f. FIG[S1 

"diagonal" cut to the parametric space in FIG. [8] 

For all four invariants under consideration, we observe 
a stronger Wfl-dependence than in Models I and II. This 
is namely due to the extended Pati-Salam running which 
contributes with larger q -coefficients than the LR stage, 
c.f., TABLElVl Moreover, unlike in FIGs.[6]and[7l three 
out of four invariants grow with lowering vr while the 
fourth one even becomes negative for vr close to the 
MSSM scale, thus, again, leading to a very characteristic 
pattern. 



C. Model IV with a sliding U{V)r x U(X)b-l scale 



The method: Finally, in Model IV, c.f. Sect. HTCl the 
unification is exact for any value of the sliding scale vbl 
below a (constant) vr in the ballpark of roughly 3 x 10^^ 
GeV, c.f., FIG. [ini Thus, as before, the main uncer- 
tainty at this level comes from the experimental error 
in OLs(}^iz) which translates into small shifts in vr. In 
what follows, we shall again vary vbl along the constant 
«s(-^^z)-error trajectories, namely, within zblcr, corre- 
sponding to the boundaries between the yellow and white 
areas within the parametric depicted in FIG.fTOl 

The results: In the two panels of FIG. [Til the four in- 
variants of our interest are depicted as functions of vbl- 
As in the case of Model III, for each of them the solid 
line corresponds to the central-value trajectory in the 
parametric space of FIG. [TOl whereas the dashed and 
dotted curves refer to the —Icr and -f-lcr-trajectories, re- 
spectively. 
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FIG. 10: The parameter space of Model IV of Sect.|ITCl The 
contours correspond to the different quality of the as{Mz) fit 
for different choice of the L-R breaking scale vr, namely, to 
1(7 (white), 2a- (yellow) and 3a (orange) values for the range 
quoted in [STJ]. 

Due to the very special nature of the sliding scale 
in this setting, all four invariants exhibit only a very 
mild VBL dependence, with the strongest effect of the 
order of few per cent observed in the LE case. This is 
because the vbl scale enters into the soft masses only 
through the slight changes in the abelian gauge couplings, 
which, however, are overwhelmed by the colour effects 
in all the other invariants. This, however, will make it 
rather difficult to distinguish this model from the MSSM, 
namely because such a discrimination is efficient only if 
more than a single invariant differs significantly from the 
mSugra value so that the intermediate scale can be inde- 
pendently constrained from more than a single quantity. 

Finally, let us comment in brief on the case where 
the U{1)r X U{1)b-l — > U{1)y breakdown is not trig- 
gered by SU{2)r doublets like above but by, e.g., SU{2)jj 
triplets. We expect that for such models the effects 
on the invariants would be similar to those expected in 
Model IV; and certainly smaller than those observed in 
Models I-III. 



D. Squark and slepton spectra. 

In FIG.[n]we plot the shapes of the MSSM squark and 
slepton spectra obtained in mSugra and in Models I, II 
and III calculated for the SPSS benchmark point, i.e. for 
Too = 90 GeV and Afi/2 = 400 GeV. This figure is to be 
understood only as an illustrative example of the differ- 
ent spectra generated in our different models. For each 
of the cases, the horizontal lines (bottom to up) corre- 
spond to TOgc (light blue), m^ (blue), tti^^c (orange), ttij^ 
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FIG. 11: The ^sL-scale dependence of the RGE invariants in 
Model IV. For practical purposes the figure has been split into 
two panels. For each of the four invariants, the colid curves 
correspond to as{Mz) fixed at its central value while the 
dashed and dotted lines refer to the — Icr and +lcr trajectories 
corresponding to roughly i)_r ~ 2 x 10^'' GeV and vr ^ ix 10^^ 
GeV, respectively; c.f. FIGfTOl 



(light orange) and rriq (purple). In order to pronounce 
the differences, the vr scale has been in all cases chosen 
very low, namely, vr ^ 10"^ GeV, and consequently vps 
in Model III is fixed to vps ^10^ GeV by gauge unifica- 
tion. The masses of the d^ and of the u'^ almost coincide 
in all the Models. Models I and II differ from the mSugra 
case namely by the smaller splittings observable in the 
squark as well as in the slepton masses, which is more 
pronounced for the latter model. However, the spectrum 
of Model III is strongly compressed due to an extended 
Pati-Salam stage which makes it rather outstanding. 

We decided not to overpopulate the figure by display- 
ing the gaugino masses which, indeed, are obtained by 
a simple rescaling Q; the SUSY-to-GUT-scale ratios of 
the relevant a's can be inferred from the evolution of the 
gauge couplings, c.f. FIGsIH FIGsH FIGs[2]and FIGslH 
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FIG. 12: The MSSM squark and slepton spectra mSugra and 
Models I, II and III calculated for the SPSS benchmark point, 
i.e. for mo = 90 GeV and M1/2 = 400 GeV. In all cases, 
VR = 10^ GeV and vps - lO'^ GeV in Model III . From 
botton to up the horizontal lines correspond to: mgc (light 
blue), TTij- (blue), mac (orange), mj^ (light orange) and rriq 
(purple). We do not show the results for model IV in this 
figure, since they are very similar to the niSugra case. 

V. DISCUSSION AND OUTLOOK 

We have studied the leading-log RGE evolution of the 
MSSM soft SUSY breaking parameters for four differ- 
ent GUT models with mSugra boundary conditions. Al- 
though all the settings are based on the unified 5*0(10) 
gauge group, they differ at the level of intermediate 
scale symmetry groups and/or particle content below 
the GUT scale. Two of the models discussed (Mod- 
els I and II of Sects. IIIBlal and IIIBlb|) . which dif- 
fer only in their beyond-MSSM field contents, feature 
an intermediate left-right symmetry which, at the level 
of precision used in this calculation, can be broken to 
SU{3)c X SU{2)l X U{1)y of the MSSM almost any- 
where between AIq and the soft SUSY-breaking scale. 
In Model III (c.f., Sect. IIIB2|) the sliding nature of the 
5t/(2)ij-breaking scale relies on an additional interme- 
diate Pati-Salam symmetry. Finally, in Model IV (see 
Sect. IIICp , the left- right symmetry is broken at a rel- 
atively high scale, but there is instead a sliding scale 
corresponding to the breaking of its U{1)b. x U{1)b-l 
remnant. All models we consider are able to accommo- 
date the neutrino data by either inverse or linear seesaw. 

The extra gauge groups and/or beyond MSSM fields 
change the evolution of the soft parameters with respect 
to the basic mSugra expectation. The invariant mass 
combinations we considered are especially suited to un- 
cover the effects of beyond-mSugra physics on the SUSY 
spectra. Remarkably, while invariants contain only a log- 



arithmic dependence on the new physics scales, their be- 
havior is qualitatively different in different models. 

In our Models I and II, the invariants LE and QU (c.f.. 
Sect, imp are always lower than the mSugra limit, while 
DL and QE are always larger. The former is a direct con- 
sequence of the LR symmetry, while the latter reflects 
mainly the shift in a{MQ) the models exhibit with re- 
spect to the MSSM expectation. Moreover, in spite of 
only a mild difference in the particle content, the invari- 
ants differ quantitatively between Model I and Model II. 

In contrast to that, in the Pati-Salem based Model III, 
LE and QU are always larger than in mSugra, with a 
rather strong dependence on the vr scale, namely due 
to the higher dimensionality of the relevant multiplets at 
the Pati-Salam stage. At the same time, in Model III, 
DL is always below the mSugra limit, while QE hardly 
varies at all as a function of vr. Finally, Model IV is an 
example of how a new scale can be effectively "hidden" 
from the RGE invariants in special constructions: De- 
spite containing a new scale potentially as low as 0(1) 
TeV, all invariants are always very close to the mSugra 
limit in this model. Technically, this is achieved by main- 
taining the beta coefficients for the SU{2)l and SU{3)c 
factors as in the MSSM all the way up to a scale close 
to AIq, while the sliding feature of the U{1)rxU{1)b-l 
breaking scale "shields" all invariants from the effects of 
the new group, with the exception of LE, which, however, 
changes only very weakly. 

It is especially interesting to compare our results 
with those obtained for minimal seesaw models within 
mSugra. Invariants for seesaw have been studied previ- 
ously in 281 |47l |49[. Type-I seesaw adds only singlets 



to the MSSM and thus, just like our Model IV, can not 
be distinguished from the pure mSugra case by means of 
the invariants only. Type-II and type-Ill seesaw, on the 
other hand, change the 5-coefficients with respect to the 
MSSM, but do not extend the gauge group. As a result, 
for minimal seesaws all four invariants are larger than 
their mSugra limit if the seesaw scale is below the GUT 
scale, as indicated by neutrino data. Thus, the invariants 
should allow to distinguish our 5O(10)-based Models I to 
III from type-II and type-Ill seesaw. 

The RGE invariants are, therefore, good model dis- 
criminators, at least in principle. However, any attempt 
to quantitatively determine the scale of a new physics 
within a particular scenario must inevitably address the 
accuracy of their calculation. Different types of errors 
need to be considered here. First, there are the errors 
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from uncertainties in the values of the input parameters. 
The largest error currently stems from the completely 
unknown msuSY, see FIG. [5] and Eq. ^. Once SUSY 
masses, indeed, have been measured, this will become ir- 
relevant and the largest error will, most likely, be lS.{as)- 

Next, the RGE invariants considered here are calcu- 
lated to the leading-log precision only. However, in some 
cases, important higher order effects such as genuine 2- 
loop corrections and 1-loop thresholds can emerge; for 
the seesaw, this was studied recently in 28|, |49[. Both, 2- 
loop running and 1-loop SUSY thresholds can, of course, 
be taken into account, but the calculation of the invari- 
ants at this level can not be done analytically. Instead, 
it requires numerical tools such as, e.g., SPheno J58ll59| 
and SARAH """" """" 



sB-lfii] 



Probably more important than the above theoretical 
considerations, eventually, will be the fact that the invari- 
ants are not directly measurable quantities. Conversion 
of the invariants into the measured sparticle masses (or 
extraction of relevant soft parameters from sparticle mea- 
surements) requires additional experimental input. In 
case of the first two generations of sfcrmions this requires 
at least a reliable measurement of tan /3 for the determi- 
nation of the D-terms. In addition, at variance with the 
situation in the minimal seecsaw models, the breaking of 
the extra gauge symmetries can potentially produce new 
D-tcrms not present in the SU{i)c x SU{2)l x U{1)y 
case. Usually it is assumed that any beyond- SM gauge 
group is broken in a "D-fiat" manner in order to avoid 
problems with tachyonic sfermions. However, since we 
wish to extract information from sfermion masses them- 
selves, it will certainly be prudent to do a combined fit 
on the new parameters instead of simply assuming D- 
flatness. 

The prospects of measuring sparticle masses at the 
LHC and, possibly, at the future ILC have been studied 
by many authors, for a detailed review sec, for exam- 
ple [6|. In general, one expects that slepton and gaugino 
masses, if within the kinematical reach of the ILC, can be 
measured at the per mill level or even better. Coloured 
sparticles, however, might be too heavy to be produced 
at the ILC. At the LHC, the precision with which spar- 
ticle masses can be measured depends strongly not only 
on the absolute scale of the SUSY masses, but also de- 
cisively on the mass ordering of the sparticles. If long 
decay chains such as q -^ x§g with x§ ^- // ^ '^'^Xi ^'^^ 
available, many SUSY masses can be measured with ac- 
curacies down to (few) percent. From the detailed studies 



of [6[ , the authors of [50| concluded that the precision of 
ILC-I-LHC combined would make it possible to see indi- 
cations for a seesaw of either type II or type HI for nearly 
all relevant seesaw scales. In an LHC-only analysis, the 
seesaw scale must be below 10^"' GeV even in favourable 
circumstances 150| or might not leave a trace in the LHC 
data at all. 

Comparing roughly the changes in spectra induced in 
the seesaw models studied in |50| with the changes ex- 
pected in our 50(10) models, we expect that a detailed, 
numerical calculation should be able to probe most, if 
not all the interesting parameter space of our models, if 
SUSY is found at the LHC and precise mass measure- 
ments are done with the help of an ILC. 

Finally, we would like to mention that the models we 
have studied in this paper have potentially also a rich 
phenomenology beyond the MSSM apart from the invari- 
ants. There are the new gauge bosons, additional Hig- 
gses, additional gauginos/higgsinos, large lepton flavour 
violation and many other effects worth studying. We 
plan to return to these questions in a future publication. 
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Appendix A: One-loop running with U(l) 



mixing 



In this appendix we give some technical details of the 
one-loop evolution of gauge couplings and soft-SUSY- 
breaking terms in Model IV of Sect. Ill CI in which ex- 
tra kinetic mixing effects, generally present in theories 
with multiple U{1) gauge factors, emerge. This, in the 
approach advocated in, e.g., J52l |. amounts to extending 
the notion of the individual gauge couplings and gaugino 
masses associated to different U{\) gauge factors to ma- 
trix forms, which, subsequently, complicates the relevant 
generalized evolution equations. 
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Gauge couplings 



b. Soft SUSY-breaking terms 



To deal with the effects of the kinetic mixing in cases 
with more than a single abelian gauge factor like in Model 
IV of Sect. Ill Cj it is convenient to work with a matrix of 
gauge couplings rather than with each of them individ- 
ually, which would require an extra RGE for the kinetic 
mixing parameters, c.f., [S^- In the C/(1)_r x U{1)b-l 
case this amounts to defining 



Neglecting for simplicity the Yukawa couplings and 
also the "Trace-terms" (denoted by S in [5l|), which in 
mSugra yield only sub-leading correction to the leading- 
log approximation used in this work, one can write the 
generalized evolution equation including the effects of the 
C/(l) mixing [52] as 



G = 



9rr 9rx 
9xR gxx 



(Al) 



where X is a shorthand notation for the canonically nor- 
malized B — L. The evolution equation can be then writ- 
ten as 



dt 



•1: 



(A2) 



where A'^ = 47r(G'G^)"^ and t = ^log(^/^o)- Here 
we have defined the relevant matrix of anomalous dimen- 
sions by 



7 



E^/Q/' 



(A3) 



where the summation is taken over all the chiral super- 
fields / in the model and Qf denotes a column vector of 
U{1)e and U{1)b-l charges of each /. 

The matching condition between such high-energy 
gauge couplings (corresponding to U{1)b. ® U{1)b-l in 
the case of our interest) and the effective-theory one (i.e., 
U{1)y oithc MSSM) at scale to then reads 



aY\to)^p^A ^{to)py 



(A4) 



where py — (y |i v |) ^^^ the coefficients of the hy- 
percharge Y in the space of the R- and B — L-charges, 

I AT. Thus, one has 



namely, Y = JfT| 



-2 ,' \-2 

9y — (9RR9XX - 9RX9XR) 



{9xx + 9xr) (A5) 



2 2 

+ 5 {Srr + 9rx) - -^^{9RR9XR + 9Rx9xx) 



Ato2 ^ --QTQMM'^G'^Qf 



dt 



(A6) 



where G is the matrix of gauge couplings, M is the gaug- 
ino mass matrix and t ~ -^ log /i//io. This is to be solved 
together with the gauge coupling (|A2p and gaugino evo- 
lution equations. The latter reads at one loop 



dt 



-M 



8tt 



{MG'^jG + G^jGM) 



(A7) 



The simphcity of the system ((A2|) . (|A6j) and (|A7)) and, 
in particular, the flavour-diagonal mSugra initial condi- 
tion, admits to write the general solution in a closed and 
compact form 



A-\t) = A-Hto)-lit-to). (A8) 
{G-'^MG-'m = i-a^'Afi/a, (A9) 

and, in particular, 



mj{t) - mjito) = 2il/2/2ap2 x (AlO) 

QJA-' [7-1 - A-'Aoj-'AoA-'] A'^Qf , 



where Aq = A{to) and A = A(t). 
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